Tn-seq of the Candida glabrata reference strain CBS138 reveals epigenetic plasticity, structural variation, and intrinsic mechanisms of resistance to micafungin.

C. glabrata is an opportunistic pathogen that can resist common antifungals and rapidly acquire multidrug resistance. A large amount of genetic variation exists between isolates, which complicates generalizations. Portable Tn-seq methods can efficiently provide genome-wide information on strain differences and genetic mechanisms. Using the Hermes transposon, the CBS138 reference strain and a commonly studied derivative termed 2001 were subjected to Tn-seq in control conditions and after exposure to varying doses of the clinical antifungal micafungin. The approach revealed large differences between these strains, including a 131 kb tandem duplication and a variety of fitness differences. Additionally, both strains exhibited up to 1000-fold increased transposon accessibility in subtelomeric regions relative to the BG2 strain, indicative of open subtelomeric chromatin in these isolates and large epigenetic variation within the species. Unexpectedly, the Pdr1 transcription factor conferred resistance to micafungin through targets other than CDR1 . Other micafungin resistance pathways were also revealed including mannosyltransferase activity and biosynthesis of the lipid precursor sphingosine, the drugging of which by SDZ 90-215 or myriocin enhanced the potency of micafungin in vitro . These findings provide insights into complexity of the C. glabrata species as well as strategies for improving antifungal efficacy. Summary: Candida glabrata is an emerging pathogen with large genetic diversity and genome plasticity. The type strain CBS138 and a laboratory derivative were mutagenized with the Hermes transposon and profiled using Tn-seq. Numerous genes that regulate innate and acquired resistance to an important clinical antifungal were uncovered, including a pleiotropic drug resistance gene (PDR1) and a duplication of part of one chromosome. Compounds that target PDR1 and other genes may augment the potency of existing antifungals.

Sustained ERK signaling promotes G2 cell cycle exit and primes cells for whole-genome duplication.

Whole-genome duplication (WGD) is a frequent event in cancer evolution that fuels chromosomal instability. WGD can result from mitotic errors or endoreduplication, yet the molecular mechanisms that drive WGD remain unclear. Here, we use live single-cell analysis to characterize cell-cycle dynamics upon aberrant Ras-ERK signaling. We find that sustained ERK signaling in human cells leads to reactivation of the APC/C in G2, resulting in tetraploid G0-like cells that are primed for WGD. This process is independent of DNA damage or p53 but dependent on p21. Transcriptomics analysis and live-cell imaging showed that constitutive ERK activity promotes p21 expression, which is necessary and sufficient to inhibit CDK activity and which prematurely activates the anaphase-promoting complex (APC/C). Finally, either loss of p53 or reduced ERK signaling allowed for endoreduplication, completing a WGD event. Thus, sustained ERK signaling-induced G2 cell cycle exit represents an alternative path to WGD.

Uncalled4 improves nanopore DNA and RNA modification detection via fast and accurate signal alignment.

Nanopore signal analysis enables detection of nucleotide modifications from native DNA and RNA sequencing, providing both accurate genetic/transcriptomic and epigenetic information without additional library preparation. Presently, only a limited set of modifications can be directly basecalled (e.g. 5-methylcytosine), while most others require exploratory methods that often begin with alignment of nanopore signal to a nucleotide reference. We present Uncalled4, a toolkit for nanopore signal alignment, analysis, and visualization. Uncalled4 features an efficient banded signal alignment algorithm, BAM signal alignment file format, statistics for comparing signal alignment methods, and a reproducible de novo training method for k-mer-based pore models, revealing potential errors in ONT's state-of-the-art DNA model. We apply Uncalled4 to RNA 6-methyladenine (m6A) detection in seven human cell lines, identifying 26% more modifications than Nanopolish using m6Anet, including in several genes where m6A has known implications in cancer. Uncalled4 is available open-source at github.com/skovaka/uncalled4.

A genome sequence for the threatened whitebark pine.

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

LARP1 senses free ribosomes to coordinate supply and demand of ribosomal proteins.

Terminal oligopyrimidine motif-containing mRNAs (TOPs) encode all ribosomal proteins in mammals and are regulated to tune ribosome synthesis to cell state. Previous studies implicate LARP1 in 40S- or 80S-ribosome complexes that repress and stabilize TOPs. However, a mechanistic understanding of how LARP1 and TOPs interact with these complexes to coordinate TOP outcomes is lacking. Here, we show that LARP1 senses the cellular supply of ribosomes by directly binding non-translating ribosomal subunits. Cryo-EM structures reveal a previously uncharacterized domain of LARP1 bound to and occluding the 40S mRNA channel. Free cytosolic ribosomes induce sequestration of TOPs in repressed 80S-LARP1-TOP complexes independent of alterations in mTOR signaling. Together, this work demonstrates a general ribosome-sensing function of LARP1 that allows it to tune ribosome protein synthesis to cellular demand.

A Genome Sequence for the Threatened Whitebark Pine.

Whitebark pine (WBP, Pinus albicaulis ) is a white pine of subalpine regions in western contiguous US and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola ) and additional threats from mountain pine beetle ( Dendroctonus ponderosae ), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short-reads of haploid megametophyte tissue and Oxford Nanopore long-reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gbp of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gbp). Approximately 87.2% (24.0 Gbp) of total sequence was placed on the twelve WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich-repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the three subclasses of NLRs (TNL, CNL, RNL). Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo assembled transcriptomes.

Comprehensive Screening of a Light-Inducible Split Cre Recombinase with Domain Insertion Profiling.

Splitting proteins with light- or chemically inducible dimers provides a mechanism for post-translational control of protein function. However, current methods for engineering stimulus-responsive split proteins often require significant protein engineering expertise and the laborious screening of individual constructs. To address this challenge, we use a pooled library approach that enables rapid generation and screening of nearly all possible split protein constructs in parallel, where results can be read out by using sequencing. We perform our method on Cre recombinase with optogenetic dimers as a proof of concept, resulting in comprehensive data on the split sites throughout the protein. To improve the accuracy in predicting split protein behavior, we develop a Bayesian computational approach to contextualize errors inherent to experimental procedures. Overall, our method provides a streamlined approach for achieving inducible post-translational control of a protein of interest.

FrameD: framework for DNA-based data storage design, verification, and validation.

MOTIVATION: DNA-based data storage is a quickly growing field that hopes to harness the massive theoretical information density of DNA molecules to produce a competitive next-generation storage medium suitable for archival data. In recent years, many DNA-based storage system designs have been proposed. Given that no common infrastructure exists for simulating these storage systems, comparing many different designs along with many different error models is increasingly difficult. To address this challenge, we introduce FrameD, a simulation infrastructure for DNA storage systems that leverages the underlying modularity of DNA storage system designs to provide a framework to express different designs while being able to reuse common components. RESULTS: We demonstrate the utility of FrameD and the need for a common simulation platform using a case study. Our case study compares designs that utilize strand copies differently, some that align strand copies using multiple sequence alignment algorithms and others that do not. We found that the choice to include multiple sequence alignment in the pipeline is dependent on the error rate and the type of errors being injected and is not always beneficial. In addition to supporting a wide range of designs, FrameD provides the user with transparent parallelism to deal with a large number of reads from sequencing and the need for many fault injection iterations. We believe that FrameD fills a void in the tools publicly available to the DNA storage community by providing a modular and extensible framework with support for massive parallelism. As a result, it will help accelerate the design process of future DNA-based storage systems. AVAILABILITY AND IMPLEMENTATION: The source code for FrameD along with the data generated during the demonstration of FrameD is available in a public Github repository at https://github.com/dna-storage/framed, (https://dx.doi.org/10.5281/zenodo.7757762).

Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but have not been considered as a feasible replacement for population-scale projects, being a combination of too expensive, not scalable enough or too error-prone. Here we develop an efficient and scalable wet lab and computational protocol, Napu, for Oxford Nanopore Technologies long-read sequencing that seeks to address those limitations. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the National Institutes of Health Center for Alzheimer's and Related Dementias. Using a single PromethION flow cell, we can detect single nucleotide polymorphisms with F1-score comparable to Illumina short-read sequencing. Small indel calling remains difficult within homopolymers and tandem repeats, but achieves good concordance to Illumina indel calls elsewhere. Further, we can discover structural variants with F1-score on par with state-of-the-art de novo assembly methods. Our protocol phases small and structural variants at megabase scales and produces highly accurate, haplotype-specific methylation calls.

Comparison of red raspberry and wild strawberry fruits reveals mechanisms of fruit type specification.

Belonging to Rosaceae, red raspberry (Rubus idaeus) and wild strawberry (Fragaria vesca) are closely related species with distinct fruit types. While the numerous ovaries become the juicy drupelet fruits in raspberry, their strawberry counterparts become dry and tasteless achenes. In contrast, while the strawberry receptacle, the stem tip, enlarges to become a red fruit, the raspberry receptacle shrinks and dries. The distinct fruit-forming ability of homologous organs in these 2 species allows us to investigate fruit type determination. We assembled and annotated the genome of red raspberry (R. idaeus) and characterized its fruit development morphologically and physiologically. Subsequently, transcriptomes of dissected and staged raspberry fruit tissues were compared to those of strawberry from a prior study. Class B MADS box gene expression was negatively associated with fruit-forming ability, which suggested a conserved inhibitory role of class B heterodimers, PISTILLATA/TM6 or PISTILLATA/APETALA3, for fruit formation. Additionally, the inability of strawberry ovaries to develop into fruit flesh was associated with highly expressed lignification genes and extensive lignification of the ovary pericarp. Finally, coexpressed gene clusters preferentially expressed in the dry strawberry achenes were enriched in "cell wall biosynthesis" and "ABA signaling," while coexpressed clusters preferentially expressed in the fleshy raspberry drupelets were enriched in "protein translation." Our work provides extensive genomic resources as well as several potential mechanisms underlying fruit type specification. These findings provide the framework for understanding the evolution of different fruit types, a defining feature of angiosperms.

Comprehensive screening of a light-inducible split Cre recombinase with domain insertion profiling.

Splitting proteins with light- or chemically-inducible dimers provides a mechanism for post-translational control of protein function. However, current methods for engineering stimulus-responsive split proteins often require significant protein engineering expertise and laborious screening of individual constructs. To address this challenge, we use a pooled library approach that enables rapid generation and screening of nearly all possible split protein constructs in parallel, where results can be read out using sequencing. We perform our method on Cre recombinase with optogenetic dimers as a proof of concept, resulting in comprehensive data on split sites throughout the protein. To improve accuracy in predicting split protein behavior, we develop a Bayesian computational approach to contextualize errors inherent to experimental procedures. Overall, our method provides a streamlined approach for achieving inducible post-translational control of a protein of interest.

Beyond assembly: the increasing flexibility of single-molecule sequencing technology.

The maturation of high-throughput short-read sequencing technology over the past two decades has shaped the way genomes are studied. Recently, single-molecule, long-read sequencing has emerged as an essential tool in deciphering genome structure and function, including filling gaps in the human reference genome, measuring the epigenome and characterizing splicing variants in the transcriptome. With recent technological developments, these single-molecule technologies have moved beyond genome assembly and are being used in a variety of ways, including to selectively sequence specific loci with long reads, measure chromatin state and protein-DNA binding in order to investigate the dynamics of gene regulation, and rapidly determine copy number variation. These increasingly flexible uses of single-molecule technologies highlight a young and fast-moving part of the field that is leading to a more accessible era of nucleic acid sequencing.

Genomic insights into metabolic flux in hummingbirds.

Hummingbirds are very well adapted to sustain efficient and rapid metabolic shifts. They oxidize ingested nectar to directly fuel flight when foraging but have to switch to oxidizing stored lipids derived from ingested sugars during the night or long-distance migratory flights. Understanding how this organism moderates energy turnover is hampered by a lack of information regarding how relevant enzymes differ in sequence, expression, and regulation. To explore these questions, we generated a chromosome-scale genome assembly of the ruby-throated hummingbird (A. colubris) using a combination of long- and short-read sequencing, scaffolding it using existing assemblies. We then used hybrid long- and short-read RNA sequencing of liver and muscle tissue in fasted and fed metabolic states for a comprehensive transcriptome assembly and annotation. Our genomic and transcriptomic data found positive selection of key metabolic genes in nectivorous avian species and deletion of critical genes (SLC2A4, GCK) involved in glucostasis in other vertebrates. We found expression of a fructose-specific version of SLC2A5 putatively in place of insulin-sensitive SLC2A5, with predicted protein models suggesting affinity for both fructose and glucose. Alternative isoforms may even act to sequester fructose to preclude limitations from transport in metabolism. Finally, we identified differentially expressed genes from fasted and fed hummingbirds, suggesting key pathways for the rapid metabolic switch hummingbirds undergo.

Nanopore sequencing identifies a higher frequency and expanded spectrum of mitochondrial DNA deletion mutations in human aging.

Mitochondrial DNA (mtDNA) deletion mutations cause many human diseases and are linked to age-induced mitochondrial dysfunction. Mapping the mutation spectrum and quantifying mtDNA deletion mutation frequency is challenging with next-generation sequencing methods. We hypothesized that long-read sequencing of human mtDNA across the lifespan would detect a broader spectrum of mtDNA rearrangements and provide a more accurate measurement of their frequency. We employed nanopore Cas9-targeted sequencing (nCATS) to map and quantitate mtDNA deletion mutations and develop analyses that are fit-for-purpose. We analyzed total DNA from vastus lateralis muscle in 15 males ranging from 20 to 81 years of age and substantia nigra from three 20-year-old and three 79-year-old men. We found that mtDNA deletion mutations detected by nCATS increased exponentially with age and mapped to a wider region of the mitochondrial genome than previously reported. Using simulated data, we observed that large deletions are often reported as chimeric alignments. To address this, we developed two algorithms for deletion identification which yield consistent deletion mapping and identify both previously reported and novel mtDNA deletion breakpoints. The identified mtDNA deletion frequency measured by nCATS correlates strongly with chronological age and predicts the deletion frequency as measured by digital PCR approaches. In substantia nigra, we observed a similar frequency of age-related mtDNA deletions to those observed in muscle samples, but noted a distinct spectrum of deletion breakpoints. NCATS-mtDNA sequencing allows the identification of mtDNA deletions on a single-molecule level, characterizing the strong relationship between mtDNA deletion frequency and chronological aging.

Multicentre genetic diversity study of carbapenem-resistant Enterobacterales: predominance of untypeable pUVA-like blaKPC bearing plasmids.

Objectives: Carbapenem-resistant Enterobacterales (CRE) are an urgent public health threat. A better understanding of the molecular epidemiology and transmission dynamics of CRE is necessary to limit their dissemination within healthcare settings. We sought to investigate the mechanisms of resistance and spread of CRE within multiple hospitals in Maryland. Methods: From 2016 to 2018, all CRE were collected from any specimen source from The Johns Hopkins Medical Institutions. The isolates were further characterized using both phenotypic and genotypic approaches, including short- and/or long-read WGS. Results: From 2016 to 2018, 302 of 40 908 (0.7%) unique Enterobacterales isolates were identified as CRE. Of CRE, 142 (47%) were carbapenemase-producing CRE with KPC (80.3%) predominating among various genera. Significant genetic diversity was identified among all CRE with high-risk clones serving as major drivers of clonal clusters. Further, we found the predominance of pUVA-like plasmids, with a subset harbouring resistance genes to environmental cleaning agents, involved in intergenus dissemination of blaKPC genes. Conclusions: Our findings provide valuable data to understand the transmission dynamics of all CRE within the greater Maryland region. These data can help guide targeted interventions to limit CRE transmission in healthcare facilities.

IVT generation of guideRNAs for Cas9-enrichment Nanopore Sequencing.

Generating high-coverage sequencing coverage at select genomic loci has extensive applications in both research science and genetic medicine. Long-read sequencing technologies (e.g. nanopore sequencing) have expanded our ability to generate sequencing data in regions (e.g. repetitive elements) that are difficult to interrogate with short-read sequencing methods. In work presented here, we expand on our previous work using CRISPR/Cas9 for targeted nanopore sequencing by using in vitro transcribed guideRNAs, with 1100 guideRNAs in a single experiment. This approach decreases the cost per guideRNA, increases the number of guideRNAs that can be multiplexed in a single experiment, and provides a way to rapidly screen numerous guideRNAs for cutting efficiency. We apply this strategy in multiple patient-derived pancreatic cancer cell lines, demonstrating its ability to unveil structural variation in "deletion hotspots" around the tumor suppressor genes p16 (CDKN2A), and SMAD4.

Similar evolutionary trajectories in an environmental Cryptococcus neoformans isolate after human and murine infection.

A pet cockatoo was the suspected source of Cryptococcus neoformans recovered from an immunocompromised patient with cryptococcosis based on molecular analyses available in 2000. Here, we report whole genome sequence analysis of the clinical and cockatoo strains. Both are closely related MATα strains belonging to the VNII lineage, confirming that the human infection likely originated from pet bird exposure. The two strains differ by 61 single nucleotide polymorphisms, including eight nonsynonymous changes involving seven genes. To ascertain whether changes in these genes are selected for during mammalian infection, we passaged the cockatoo strain in mice. Remarkably, isolates obtained from mouse tissue possess a frameshift mutation in one of the seven genes altered in the human sample (LQVO5_000317), a gene predicted to encode an SWI-SNF chromatin-remodeling complex protein. In addition, both cockatoo and patient strains as well as mouse-passaged isolates obtained from brain tissue had a premature stop codon in a homologue of ZFC3 (LQVO5_004463), a predicted single-zinc finger containing protein, which is associated with larger capsules when deleted and reverted to a full-length protein in the mouse-passaged isolates obtained from lung tissue. The patient strain and mouse-passaged isolates show variability in virulence factors, with differences in capsule size, melanization, rates of nonlytic expulsion from macrophages, and amoeba predation resistance. Our results establish that environmental strains undergo genomic and phenotypic changes during mammalian passage, suggesting that animal virulence can be a mechanism for genetic change and that the genomes of clinical isolates may provide a readout of mutations acquired during infection.

Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but to date have not been considered as feasible replacement at scale due to a combination of being too expensive, not scalable enough, or too error-prone. Here, we develop an efficient and scalable wet lab and computational protocol for Oxford Nanopore Technologies (ONT) long-read sequencing that seeks to provide a genuine alternative to short-reads for large-scale genomics projects. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the NIH Center for Alzheimer's and Related Dementias (CARD). Using a single PromethION flow cell, we can detect SNPs with F1-score better than Illumina short-read sequencing. Small indel calling remains to be difficult inside homopolymers and tandem repeats, but is comparable to Illumina calls elsewhere. Further, we can discover structural variants with F1-score comparable to state-of the-art methods involving Pacific Biosciences HiFi sequencing and trio information (but at a lower cost and greater throughput). Using ONT based phasing, we can then combine and phase small and structural variants at megabase scales. Our protocol also produces highly accurate, haplotype-specific methylation calls. Overall, this makes large-scale long-read sequencing projects feasible; the protocol is currently being used to sequence thousands of brain-based genomes as a part of the NIH CARD initiative. We provide the protocol and software as open-source integrated pipelines for generating phased variant calls and assemblies.

Treatment of Sindbis Virus-Infected Neurons with Antibody to E2 Alters Synthesis of Complete and nsP1-Expressing Defective Viral RNAs.

Alphaviruses are positive-sense RNA viruses that are important causes of viral encephalomyelitis. Sindbis virus (SINV), the prototype alphavirus, preferentially infects neurons in mice and is a model system for studying mechanisms of viral clearance from the nervous system. Antibody specific to the SINV E2 glycoprotein plays an important role in SINV clearance, and this effect is reproduced in cultures of infected mature neurons. To determine how anti-E2 antibody affects SINV RNA synthesis, Oxford Nanopore Technologies direct long-read RNA sequencing was used to sequence viral RNAs following antibody treatment of infected neurons. Differentiated AP-7 rat olfactory neuronal cells, an in vitro model for mature neurons, were infected with SINV and treated with anti-E2 antibody. Whole-cell RNA lysates were collected for sequencing of poly(A)-selected RNA 24, 48, and 72 h after infection. Three primary species of viral RNA were produced: genomic, subgenomic, and defective viral genomes (DVGs) encoding the RNA capping protein nsP1. Antibody treatment resulted in overall lower production of SINV RNA, decreased synthesis of subgenomic RNA relative to genomic RNA, and suppressed production of the nsP1 DVG. The nsP1 DVG was packaged into virus particles and could be translated. Because antibody-treated cells released a higher proportion of virions with noncapped genomes and transient transfection to express the nsP1 DVG improved viral RNA capping in antibody-treated cells, we postulate that one mechanism by which antibody inhibits SINV replication in neurons is to suppress DVG synthesis and thus decrease production of infectious virions containing capped genomes. IMPORTANCE Alphaviruses are important causes of viral encephalomyelitis without approved treatments or vaccines. Antibody to the Sindbis virus (SINV) E2 glycoprotein is required for immune-mediated noncytolytic virus clearance from neurons. We used direct RNA nanopore sequencing to evaluate how anti-E2 antibody affects SINV replication at the RNA level. Antibody altered the viral RNAs produced by decreasing the proportion of subgenomic relative to genomic RNA and suppressing production of a previously unrecognized defective viral genome (DVG) encoding nsP1, the viral RNA capping enzyme. Antibody-treated neurons released a lower proportion of SINV particles with capped genomes necessary for translation and infection. Decreased nsP1 DVG production in antibody-treated neurons led to lower expression of nsP1 protein, decreased genome capping efficiency, and release of fewer infectious virus particles. Capping was increased with exogenous expression of the nsP1 DVG. These studies identify a novel alphavirus DVG function and new mechanism for antibody-mediated control of virus replication.

Long read mitochondrial genome sequencing using Cas9-guided adaptor ligation.

The mitochondrial genome (mtDNA) is an important source of disease-causing genetic variability, but existing sequencing methods limit understanding, precluding phased measurement of mutations and clear detection of large sporadic deletions. We adapted a method for amplification-free sequence enrichment using Cas9 cleavage to obtain full length nanopore reads of mtDNA. We then utilized the long reads to phase mutations in a patient with an mtDNA-linked syndrome and demonstrated that this method can map age-induced mtDNA deletions. We believe this method will offer deeper insight into our understanding of mtDNA variation.